DOES THE GLUON SPIN CONTRIBUTE IN A 
GAUGE-INVARIANT WAY TO NUCLEON SPIN? 



Pervez Hoodbhoy* and Xiangdong Ji 
Department of Physics 
University of Maryland 
College Park, Maryland 20742 
(UMD PP#00-008 DOE/ER/40762-190 August 1999) 



Abstract 



Although the matrix element of the gluon spin operator in nucleon helicity 
states is known to be independent of some special choices of gauge, we show 
that it is not invariant under a general gauge transformation. We find that 
there exists a simple means of obtaining the matrix element in a different 
choice of gauge from a calculation made in one specific gauge. Similar con- 
clusions hold for other manifestly gauge-dependent operators present in the 
QCD angular momentum operator. 
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In quantum chromodynamics (QCD), the gluon spin operator, S g , can be denned as, 



S g = Jd 3 xExA, (1) 

where we have suppressed the color indices, E is the color electric field and A the gauge 
potential. The above operator is manifestly gauge dependent, and therefore one expects 
that its matrix element in a physical state depends on gauge choice. Indeed, in Ref. it 
was shown that the matrix element of S g in a quark helicity state is different in the axial 
(A + = 0) and covariant gauges. On the other hand, it was claimed in Ref. || that this 
matrix element in a helicity state is invariant under a gauge transformation. The goal of 
this paper is to clarify the relationship of the two results. In particular, we point out that 
the proof in Ref. covers only a special class of gauge transformations and the matrix 
element of S g is in fact gauge dependent. We then make a number of pertinent observations 
about comparing calculations in different gauges. In particular, we demonstrate that there 
exists a simple means of getting the answer in a different choice of gauge from calculations 
made in one specific gauge. 

We start with general comments about calculations in quantum gauge theories. Because 
of the gauge symmetry and the infinite number of degrees of freedom in fields, it is always 
necessary to choose a gauge for calculating Green's functions and physical matrix elements, 
perhaps with the exception of lattice gauge theory in which fields are assigned at a finite 
number of spacetime points. Indeed, in a canonical quantization of gauge theories, one has 
to choose a gauge at the very outset. In the path integral formulation, a gauge choice is 
conveniently made by selecting a set of gauge conditions, 

F a (A) = w a (x) , (2) 

where the index a runs over the number of generators of the gauge group (and will be 
omitted in the remaining equations). In addition, one can choose an arbitrary weighting 
functional G(w) to integrate over the auxiliary field w. Without loss of generality, one can 
assume J[Dw]G(w) = 1, where D denotes the functional integration measure. Hence, the 
generating functional for Green's functions in QCD reads ||, 

Z(J) — — J [D(f>] [Dw]A F (A)S(F(A) - w)G(w) exp (iS + i J <fx[J ■ A + fjip + tpr]}^j , (3) 



where <ft denotes a collection of fields A, t/>, and ip, and S is the canonical QCD action, 



(4) 



S = J d 4 x 

and the normalization constant M is 

M = J [D<p}[Dw]A F (A)5(F(A) - w)G{w)e iS , (5) 

so that Z[0] = 1. The Faddeev-Popov determinant Ap{A) is defined such that 

A F {A) J \Du]5{F{ w A) -w) = l (6) 
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where U A is a gauge transformation of A and [Duo] integrates over the gauge group at every 
point of spacetime. 

To calculate the physical matrix element (f\0(<f)\i) of an operator 0(0), one can start 
with the following Green's function, 

(0\T(J f (x)O(y)jK0))\0)\ F , G 
= ±jf [D$\[Dw]A F {A)5{F{A) - w)G(w)J f (<j ) )0( ( j ) )j}( ( f ) )e iS , (7) 

where Jij are the interpolation fields or currents for the initial and final physical states. 
According to the Lehmann-Symanzik-Zimmermann reduction formula |J, (f\0(<f))\i) is just 
the residue of the Green's function at the poles corresponding to external physical states, 
modulo the coupling constants Z { j defined as (0| J^(0)|z, /) = J~Zij (where a possible 
Lorentz structure is suppressed). Gauge invariance of the physical matrix element means 
that the residue depends on a choice of F(A) and G(w) only through the couplings Z it f 
when the interpolating currents are gauge dependent. 

Using methods available in textbooks one can show that the Green's function 
(0|T( Jf(x)0(y) j}(0))\0) \f,g is independent of the choice of F(A) and G(<j>) provided that 
Ji(4>), Jf{4>)i an d 0(4>) are gauge invariant. Indeed, consider a different choice of F and G 
such that 

Ap{A) J [Du]S{F{ u A) -w) = l; J [Dw]G(w) = 1 (8) 

Multiplying both unities to the right-hand-side of Eq. (|7|), we get 

(0\T(J f (x)O(y)jj(0))\0)\ FiG 
= 1 j [D<j>][Dw][Dw][Duj] A F {A)8{F{A)-w)G{w)J f {<j>)0{<t>)Jt^) 

x A F (A)G(w)5(F( ul A) -w)e iS 
= jfj [D<p][Dw][Dw][Du] A F {A)5(F{ U '' 1 A) -w)G{w)J f (^ 1 (j))0{ ul ~ 1 (f)) 

xj}^ 1 cj))Ap{A)G{w)5{F{A) - w)e iS 
= jfj [D<j>][Dw][Dw] G(w)J f r<j ) )0(^)j}(-y)A p (A)G(w)5(F(A)~w)e iS 

= ^J [D4>][Dw] J f {4>)0{^{^)A F {A)G{w)8{F{A) - w)e tS } (9) 

where after the second equal sign, we have made the gauge transformation <p — ^ 1 <fi and 
used the fact that the measure [D<p} and determinants A F (A) and A F (A) are invariant under 
the transformation. After the third equal sign, we have integrated out the [Du] by using the 
5-function constraint F(^ A) = w. For a fixed A and w, there is a special Uq that fulfills 
the constraint, and hence c^o is a specific function of them: u> = u~ 1 (A, w). After the fourth 
equal sign, we have used the assumption that J% t f(<f>) and 0(<p) are gauge invariant. This 
makes the functional integration over w trivial. The final line establishes gauge independence 
of the Green's function: (0\T(J f (x)O(y)j}(0))\0)\ FjG = (0\T(J f (x)O(y)j}(0))\0)\ FG . 

From the above discussion, we see that gauge invariance of a matrix element is guaranteed 
when the observable under consideration, defined as a functional of physical fields in the 



3 



path integral formalism, is invariant under any gauge transformation u, including arbitrary 
dependence on the gauge potential A itself. The notion of A-dependent gauge transformation 
can be found in textbooks, but it is usually discussed in the context of classical gauge 
theory. In a canonically-quantized theory, an A-dependent to can no longer be considered 
as a transformation "parameter" , but rather a quantum operator in Hilbert space. Because 
the commutation relations between u and the fundamental fields depend on the choice of 
transformation itself, it is difficult to consider the most general gauge transformations in 
canonically-quantized theories. In fact, as we have alluded to before, canonical quantization 
of a gauge theory can only be carried out after a choice of gauge, and different gauge 
choices often lead to different quantum Hilbert spaces. When comparing calculations in 
different gauges, one usually compares the final results only— not the actual physical states 
and operators in those gauges. In particular, we know of no discussion in the literature 
about how to transform a quantum operator from one gauge to another. 

Fortunately, in the path-integral formulation, all field variables are treated as classical, 
and a transformation from one gauge to another, performed with an A-dependent u, can 
be easily handled by a change of integration measure and the associated Jacobian (the 
Faddeev- Popov determinant). A decisive test of gauge invariance of an observable can be 
done with a general gauge transformation: an 00(A) solving F( W A) = w for any A satisfying 
F(A) = w. It can happen that certain matrix elements are invariant under a special class 
of gauge transformations; in particular, those that keep the Faddeev-Popov determinant 
invariant. However, if they cannot pass the general test above, one cannot claim their total 
gauge independence. 

In Ref. 0, Chen and Wang claim to have shown that the matrix element of S g in a 
helicity eigenstate is gauge invariant. What they have actually shown is that the matrix 
element is the same in gauges F(A) = w and F( u A) = w, where u> is an A-independent gauge 
transformation parameter. An outline of the proof goes as follows (using gauge-invariant 
interpolating operators for external states): 

(0\T(J f (x)S g (y)jJ(0))\0)\ F( , A) , G 

= ^J [D<p}[Dw}A F (A)5(F(-A)-w)G(w)J f (<f ) )SMJK<Py S 

= i J [D<j>] [Dw] A F (A)5(F(A) - w)G(w) J f {^)S g {^ 4>)Jl{4>Y S (10) 

where after the second equality one has made a gauge transformation tu -1 , and all the u 
dependence now appears in S g . Write 

^rVi^.w+^w) . (ii) 

Chen and Wang showed that 5S g (<j), u) has a zero matrix element in a helicity state and so, 

(0\T(J r (x)S g (y)j}(0))\0)\ F{ , AlG = (0\T(Ji{x)S B (y)f f (0)m\F(A),a • (12) 

As we have demonstrated above, Eq.([Tj|) is not sufficient to guarantee that that matrix 
element remains invariant under an arbitrary new gauge condition G(A) = w. Indeed, 
under A-dependent gauge transformations, the proof in Ref. no longer applies. 
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As an example of A-dependent gauge transformation, we compare calculations in the 
covariant gauge and light-cone gauge (A + = 0). The vacuum matrix element of a general 
operator O (which can be a time-ordered product of operators at several different spacetime 
points) in the covariant gauge can be written as, 

(0|6|0)|i = jfj [D<f>] [Dw)A F (A)5(d ■ A - w)0{<p) exp(-iw 2 /2X) exp(zS) , (13) 

where we have included the standard Gaussian weighting function. The same matrix element 
in light-cone gauge A + = can be expressed as, 

(0|O|0)U = jfj [ J D0]A F (A)5(A + )O(0)exp(^) . (14) 

Using the same method used in deriving Eq.(9), we can show that the matrix element in 
two different gauges is related by the following equation: 

(o|6(^)|o>u = (o|6(^)|o>u, (is) 

where u is a gauge transformation which brings a gauge configuration to the A + = gauge, 

co = -±A + . (16) 

Equation ( |T5D is important. It relates any matrix element (or Green's function) in the axial 
gauge to that of a gauge-transformed operator in the covariant gauge! For instance, using 
the relation one can easily recover the gluon or quark field propagators in the light-cone 
gauge from those in the covariant gauge. We emphasize here that similar relations between 
any two independent gauges can be derived and they are entirely nonperturbative. The 
relationship between the covariant and axial gauge calculations has also been studied recently 
by Joglekar and collaborators ||. In those studies, generalized BRST transformations are 
used to connect different gauge choices. 

As an application of Eq.(^), we consider the difference of the S g matrix elements in the 
covariant and light-cone gauges, 

(f\S g \i)\ A - (f\S g \i)\ L = (f\S g ( u A) - S g (A)\i)\ L , (17) 

which can be calculated completely in the covariant gauge. Given u in Eq.flTBD, one cannot 
show, perturbatively or nonperturbatively, that the right-hand side vanishes when the ex- 
ternal states have definite helicity. In fact the above equation can be used to directly check 
the one-loop calculation presented in Ref. [0. Consider an "on-shell" quark in the state of 
momentum p M and helicity 1/2. We notice that the quark fields are not gauge invariant. 
However, the difference between the original and the gauge-transformed quark fields has no 
perturbative quark pole; therefore, Eq. (|T7|) still applies. At one-loop order we find, 

( P+ i|^)|p + i) = CF |l„g), (18) 

where 6S g = S g ( w A) - S g (A), C F = (iV c 2 - l)/(2iV c ) with N c the number of colors, and 
Q 2 and fi 2 are the ultraviolet and infrared cutoffs, respectively. The above result is exactly 
what we found in Ref. JTJ. 
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We have also found that in the one-loop calculation the S g matrix element in the covariant 
gauge is independent of the Feynman parameter A. By studying the relations between matrix 
elements in d^A^ = w and d^A^ = w/V\ gauges, we can show the A- independence directly. 
Notice that the different choice in A, like the gauge transformations considered in Ref. [0, 
leaves the Faddeev- Popov determinant invariant. 

To summarize, we demonstrated that the general gauge-independence of a physical ma- 
trix element must be checked in the path integral formalism through gauge-field-dependent 
gauge transformations. Using this, we showed that the gluon spin contribution to the nu- 
cleon spin is indeed gauge- dependent. The conclusion also applies readily to other gauge- 
dependent operators considered in Ref. 0. We derived a general relation between matrix 
elements in the covariant and axial gauges. Using the relation, we calculated the one-loop 
difference of the gluon spin contribution to the quark spin in the two gauges. The result 
confirms the explicit calculation presented in Ref. 0. 
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